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Abstract 

We suggest the Hamiltonian approach for fluid mechanics based 
on the dynamics, formulated in terms of Lagrangian variables. The 
construction of the canonical variables of the fluid sheds a light on 
the origin of Clebsh variabes, introduced in the previous century. The 
develoted formalism permitts to relate the circulation conservation 
(Tompson theorem) with the invariance of the theory with respect to 
special diffiomorphisms and estublish also the new conservation laws. 
We discuss also the difference of the Eulerian and Lagrangian descrip- 
tion, pointing out the incompeteness of the first. The constructed 
formalism is also applicable for ideal plasma. We conclude with sev- 
eral remarks on the quantization of the fluid. 

1 Introduction 

At the present time almost all fundamental physical phenomena could be for- 
mulated in the frameworks of either classical or quantum mechanics. That 
means that these phenomena admits the Hamiltonian description, which due 
to its long history developed many powerful methods of analysis of the gen- 
eral properties of evolution of the systems and the tools for the solutions of 
partial problems. The Hamiltonian formalism also provides the unique way 
of transition from classical to quantum description of the systems. 

In this respect the fluid mechanics stands aside (in spite of its name) from 
orthodox mechanics. The reasons for that is not only the infinite number of 
degrees of freedom of fluid. We have already learned how to formulate the 
theory of classical and quantum fields. The main difference between con- 
ventional field theory and the fluid is that in the first case we can speak 
about the dynamics of the field at one point in space (which of course in- 
teracts with the field at the neighbouring points), while in the case of the 
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fluid, describing the interaction of tlie neiglibouring particles which consti- 
tute the fluid we are loosing its position in the space due to the motion of 
fluid. In the same time the objective of usual problems of hydrodynamics is 
to define the velocity, density and a thermodynamical variable (pressure or 
entropy) as the functions of the coordinates x and time t for the appropriate 
boundary conditions and/or initial data IQ. The similar problems also ap- 
pear for magnetohydrodynamics dealing with sufficiently dense plasma [@]. 
For the developing of the Hamiltonian formalism we need to start with more 
detailed description based initially on the trajectories of the particles which 
constitutes the fluid or plasma. This description is especially important for 
plasma, because the fundamental electromagnetic interaction could be for- 
mulated only in terms of the trajectories of the charges. Needless to say 
that some aspects of this approach was extensively studied in the series of 
papers by J.Marsden, A.Weinstein, P.Kupershmidt, D.Nolm , T.Ratiu and 
CLevermore especially in the context of stability problem for ideal fluid. 
In contrast to these papers here we deliberately avoid, when it is possible, the 
use of language of modern differential geometry to make our paper accessible 
to physicists. 

Of course not all properties of fluid could be formulated in the frameworks 
of the Hamiltonian approach. For example, we leave open the question of 
the energy dissipation, viscosity et cetera. 



2 The Lagrangian equations of motion 

In fluid mechanics there are two different pictures of description. The first, 
usually refereed as Eulerian, uses as the coordinates the space dependent 
fields of velocity, density and some thermodynamic variable. The second, 
Lagrangian description, uses the coordinates of the particles x{^i,t) labeled 
by the set of the parameters ^j, which could be considered as the initial po- 
sitions ^ = x{t = 0) and time t. These initial positions ^ as well, as the 
coordinates x{C,i,t) belong to some domain D C R^. In sequel we shall con- 
sider only conservative systems, where the paths of different particles do not 
cross, therefore it is clear that the functions x{^i,t) define a diffeomorphism 
oi D C and the inverse functions C,{xi,t) should also exist. 
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^^•(--^).^.(,,.) =^r (1) 
The density of the particles in space at time t is 

p{x, t)= f d3^pote)5(^ - f te, t)), (2) 



where po(0 is the initial density at time t = 0. The velocity field w as a 
function of coordinates x and t is: 

v{xi,t) = x{^{xi,t),t), (3) 

where ^(a;,t) is the inverse function (|l|). The velocity also could be written 
in the following form: 

or 



p{xi,t)v{xi,t) = J d ^po{^i)x{^i,t)S{x - x{^i,t)). (5) 
Let us calculate the time derivative of the density using its definition (H) : 

p{xi,t)= j d^ipo{ii)^^5{x-x{ii,t)) 
d 

= --^pixi,t)v{xi,t) (6) 

In such a way we verify the continuity equation of fiuid dynamics: 

p{xi, t) + d[p{xi, t)v{xi, t)) = 0. (7) 

Using the coordinates x{C,i,t) as a configurational variables we can con- 
sider the simplest motion of the fiuid described by the Lagrangian 

L^fd^^n^lMl, (8) 
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The equations of motion which follow from (j^) apparently are 

mf(ei,t)=0 (9) 

Now let us find what does this equation mean for the density and velocity 
of the fluid. For that we shall differentiate both sides of with respect to 
time 

d 



(10) 



The first term in the r.h.s. of (|l^) vanishes due to the equations of motion(P) 
and transforming the second in the same way, as we did in (10) we arrive at 



d d 

—p{Xi,t)v{Xi,t) + —(^p{Xi,t)v{Xi,t)VkiXi,t) 

Let us rewrite (0) in the following form: 

d 







(11) 



v{Xi,t) p{Xi,t) + —(^p{Xi,t)Vk{Xi,t)^ 



d 



+p{Xi,t) v{Xi,t) +VkiXi,t)- — v{xi,t) 
L dXk 



0. 



(12) 



The first term in (|T^) vanishes due to the continuity equation, while the 
second gives Euler's equation in the case of the free flow: 



d 

v{Xi, t) + Vk{Xi, t)- — v{xi, t) =0 
OXk 



(13) 



To move further we need to introduce into the Lagrangian (^ the "po- 
tential energy" term which will give rise to the internal pressure field in 
Euler's equation. As we have mentioned above, the functions x{C,i,t) define 
a diffeomorphism,in therefore the Jacobean matrix 



dXj{^i,t) 

dik 



(14) 
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is nondegenerate for any C,i and t. The integral (Q) may be expressed via the 
Jacobean determinant of A^^i, t): 



p{x, t) 



(15) 



Assume for the simphcity that the initial density po(^i) is uniform in D C E?^ 
and effectively normalize the density field p{^i,t) by putting po(^i) = 1 (one 
particle in the elementary volume). Then (|15D becomes 



detA{^i,t) 



Consider as the "potential energy" the functional of detA: 



d^^ 



f{detA{i,,t)) 



Now the equations of motion become 

d 



'{ii,t)f'{detA)detA \ = 



(16) 



(17) 



(18) 



Substituting Xj{^i,t) from (|T8D to the equation(|TOD and acting as we did in 
the derivation of the equation (pTsP, we obtain 



Vj{Xi,t) - 



_d_ 

dxj 



f\detA{i,,t)) 



= 
(19) 

It is now obvious that if we identify the —1/mf'idetAi^i, t)) ^ with 

pressure p{xi,t), the equation (^) takes the form of usual Euler equation 
without viscosity: 



p{Xi,t) l^+Vk{Xi,t)-^jVj{Xi,t) 



d 



dx-i 



■p{xi,t). 



(20) 



Note, that the pressure p{xi,t) which appeared here is not due to an ex- 
ternal force, but the result of interaction between particles. The further 
complication of the Lagrangian is not necessary for the moment. One com- 
ment nevertheless should be done. No modification of the Lagrangian will 
give us the Navier-Stokes equation because the later includes the dissipation 
effects, which is not time symmetric... 
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3 Hamiltonian formalism 



The Hamiltonian formalism of classical mechanics, which among other adven- 
tures gives us the unique way for construction of the quantum theory. In this 
respect the hydrodynamics stands aside because its fundamental variables - 
local velocity, density and thermodynamical function does not admit an im- 
mediate Hamiltonian interpretation, in spite of quite often in the text-book 
one can meet the terms like "momentum density" (see e.g. 

We shall introduce the canonical variables according to the Lagrangian 
(p!?!) , so our canonical coordinates will be the functions x{C,i, t) and its conju- 
gated momenta are defined as the derivatives of the Lagrangian with respect 
to the velocities x{^i,t): 

SL 

pte,t) = 7T--— = mf(e.,t), (21) 

The Hamiltonian is given by the Legendre transform of the Lagrangian: 

H = jd'i (^/(6,t) + /(rfetA(e„t))) . (22) 
The canonical Poisson brackets is defined by 

fete),^fc(C')} = M'te-C') (23) 



Apparently the Poisson brackets (|23|) and the Hamiltonian ( P2D define the 
equations of motion for the canonical variables, which are equivalent to 
the Lagrange equations. The phase space F of the fluid is formed by 
x{^i,t),'p{^i,t). What is missing at the moment is the space interpretation 
of the variables x{^i, ) and p{^i). To get one let us introduce the new objects 
using the same averaging, as we used in the previous section 

lix) = I d'mi)Six-xm = pix)p{^,{x)). (24) 

This new object will serve as a part of the coordinates in the phase space 
r, therefore we need to calculate its Poisson brackets, using (p3D. The result 
has the following form: 

{lj{xi),lk{yi)} = + ljiyi)-^]6{x - y). (25) 
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The Poisson brackets was introduced geometrically as "the hydrodynamic- 
type brackets" long ago in the papers |0 without physical derivation. The 
present discussion reveals the origin of these brackets. The commutation 
relation ( ]25|) coincides with algebra of 3-dimensional diffeomorphisms, where 
the lj{xi) serves as the generators. In other words, with the help of these 
generators we can realize the finite diffeomorphism Xj — > (f)j{xi) of any x- 
dependent dynamical variable. It should be mentioned that the group of 
diffeomorphisms, generated by lj{xi) is not a gauge symmetry in case of fluid 
mechanics, as it is in the case of e.g. relativistic string or membrane. In the 
same time in fluid mechanics there is an infinite dimensional symmetry with 
respect to special (i.e. volume preserving) diffeomorphisms SDiff, which 
will be considered in the next section. 

As we mentioned above, the lj{xi) is only a part of x-dependent coor- 
dinates in the phase space F, one may consider it as the "x-dependent mo- 
menta" . Now we should define the corresponding "x-dependent coordinates" . 
For example we may choose for this role the functions ^j(xj), which are the 
inverse to Xj{^i) functions Differentiating the first equation (|1|) with 
respect to x we obtain: 



= Sj, (26) 



in other words, the matrix 



is inverse to Al.{C,i{xk)). Therefore , from (jLOj) follows that 

p(xj) = deta{xi) (28) 

To simplify the calculation of the Poisson brackets we can express C,j{xi) in 
the following form: 

^^■^"^^^ /rf3^5(f-f(6,t)) ■ ^^^^ 



From (^) we easily obtain 



{U^^),Uy^)}=0 (30) 
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The calculation of the Poisson brackets between ^j{xi) and lj{xi) is more 
involved, but the result is simple: 

{lM^.Uy^)} = (31) 

In such a way we have constructed the set of x-dependent coordinates 
{lj{xi),^j{xi)) in the phase space F, but the transformation 

{Pji^i),Xj{^i)) {lj{xi),^j{xi)) (32) 

is not canonical. The set of canonical x-dependent coordinates in F could be 
obtained in the following way. Let us multiply both sides of (pTD by matrix 

Al{Uxk)m^^),Uyi)} = -StS{S-y), (33) 

where we have used (PB]). Further, due to the relation (|5DD, we can put 
Ai^{^i{xk)) inside the brackets and obtain: 

{7rUxi),^k{yi)} =Si6{x-y), (34) 

where 

T^miXi) = -Al^{^i{Xk))lj{Xi) (35) 

By direct calculation we also obtain 

{7Cm{Xi),TTk{yi)} = 0, (36) 

so the set {'n'm{xi),^k{yi)) is formed by the canonical variables. In terms of 
these canonical variables the generators of the group of diffiomorphisms lj{xi) 
has the following form: 

lj{xi) = ~~^^^'^k{xi). (37) 

For the Lagrangian considered, from ( |21| ) and (0) follows that lj{xi) is given 
by 

lj{xi) = mp{xi)vj{xi) (38) 

and the representation (^) is very similar to Clebsh parametrization |^ of 
the velocity. The distinction of ( p7| ) from the original Clebsh paramerization 
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is the appearance of three "potentials", instead of two, for the 3- dimensional 
fluid. The reason of this difference will be discussed in the end of this section. 

The above given construction of the canonical x-dependent variables is 
the most simple one. For some problems another set might be useful. Let us 
construct one denoting by Cj{xi) the canonical coordinates: 



Calculating the Poisson brackets of Cji^i) ^"^^ ^ji^i) we obtain: 

d 

{lj{xi),Ck{yi)} = (k{xi)—6{x-y). 



(39) 



(40) 



From (^) we conclude that lj{xi) should have the following representation 



d 

where i]k{xi) is the variable, canonically conjugated to Cji^i)'- 

{Vjixi)Xkiyi)} = 6jk5ix-y). 



(41) 



(42) 



The canonical Hamiltonian (^) should now be expressed in the terms of 
new variables {'n'm{xi),^k{yi))- Indeed, let us insert the unity 

1 = J d^x6{x - x{^i)) (43) 

into the integrand ( P2| ) and change the order of integration: 

H = Jd'xJ d^^ (^^"(6, t) + fidetA^^ 6{x - f (6)), (44) 



Fulfilling the integration over ^ with the help of (|]) and (24) we obtain: 



H 



d^x 



/^(x.)+p(x.)/(— -))). 



2mp{xi] 



(45) 



Making use of (|37|) we can express H via canonical variables {'Km^Xij^^kiyi))'- 
1 dik{xi) dim{xi 



H 



d'^x 



2mp{xi) dxj dxj 



-7rk{xi)nmixi) + V{p{xi)) 



(46) 



where we have introduced a notation V{p{xi)) for the "potential" part of 
the energy. This term represents the internal energy of the fluid and it 
should vanish (up to inessential constant) for homogeneous density distri- 
bution p(xj) = po- A phenomenological expression for V{p{xi)) could be 
written as follows 0: 

Vipix,)) = ^{5p{x,)f + \{yp{x,)f + ... , (47) 

where 6p{xi) is the deviation of the density from its homogeneous distribu- 
tion: 

Sp{xi) = p{xi) - Po. (48) 

The first term in (|47|) is responsible for the sound wave in the fluid (/t is the 
velocity of sound), the second term in (|47| ) describes the dispersion of the 
sound waves. 

In order to reveal the relation of the Hamiltonian flow, generated by (^) 
with geodesic flow [Q] let us introduce the metric tensor gjk{xi): 

g,,{x,) = AJ{Ux,))A^{Ux,)). (49) 

With this notation (HBI) takes the form: 



H= d'x 



fx ( -^y^/™(x,)vrfc(a;,)7r^(x,) + V{^^)] . (50) 




2"^ 'ig[x, 

The metric tensor with upper indices g^^{xi) denotes, as usually the inverse 
matrix and ^ 

g{xi) = detgjkixi) = (51) 

p^Xi) 

The representation (|50D permit us to consider the hydrodynamics as the 
geodesic flow on the dynamical manifold with metric gjk{xi) and many gen- 
eral properties of the hydrodynamics could be derived from this fact (see e.g. 

We shall complete this section with one important comment concerning 
our description of hydrodynamics. In the present discussion the phase space 
of the 3-dimensional fluid is 6- dimensional, as it naturally follows from our 
approach. This could be compared with recent paper M, where (in our 
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notations) only p{xi) and Vj{xi) are regarded as the coordinates in the phase 
space (this point of view of could be found also in different text books and 
articles) see for example ||^, [0, JTD . 



According to the conventional point of view the state of the fluid is deter- 
mined by its velocity and density and therefore all other variables like ours 
x{^i) are not needed. Indeed, solving the Euler equations of motion we can 
express the velocity v{xi,t) and the density p{xi,t) at time t via the initial 
data V(Xi,0) and p(^Xi,0) . Then, from the definition of v{xi,t) follows 

v{x{^,t),t)) = S{^,t). (52) 

(here we have suppressed the indices of the arguments for brevity). Ap- 
parently we can solve these equations with respect to x{^,t), provided we 
know the initial data x{^, 0). Therefore it may seem that the variables x{^, t) 
are unnecessary, as it could be obtained through the others. But the initial 
data x{^,0) are half of the canonical variables in Hamiltonian formalism. 
So, in our approach we indeed need more variables. These additional vari- 
ables provide the complete description of the fluid in a sense that solving the 
equations of motion we define not only the v{xi,t) and p{xi,t) but also the 
trajectories of the particles which could not be obtained in the conventional 
formalism. The situation is analogous to the rigid body rotation: here the 
phase space F is formed by 3 angles and 3 components of the angular mo- 
mentum Jj. As the Hamiltonian depends on the components of Jj only, we 
can consider separately the evolution of the angular momentum. This is in- 
complete description for this mechanical system. The complete one certainly 
should includes the evolution of the angles, which define the location of the 
rigid body in space. 

In the case of fluid dynamics the complete description is to be done in the 
6-dimensional phase space T formed by x{^) and p(^) (or ^(x) and /(x)). If, as 
it is usually the case, the Hamiltonian depends only on l{x) and p(x), partial 
description in terms of the velocities and densities considered as "relevant" 
variables is possible. This means that we do not care about the evolution 
of the coordinates of the phase space which are considered as inessential. 
The "relevant" part of the coordinates does not necessarily form a simplectic 
subspace in F. The nondegenerate Poisson brackets in F could become de- 
generate on the subset of F, corresponding to the "relevant" variables. This 
is indeed the case in the rigid body and in the conventional fluid dynamics. 
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In the first case the degeneracy of the algebra of Poisson brackets for the "rel- 
evant" variables — angular momentum is well known. Its center element is 
Casimir operator of the rotation group. In the case of fluid dynamics the al- 
gebra of the "relevant" variables — velocities and densities has the following 
form: 

{vj{xi),Vk{yi)} = -'^kVj{xi))6{x - y) 

{vj{xi),p{yi)} =Vj6{x-y) 

{pix,),p{y,)} =0 (53) 



This algebra could be obtained from equations (pS]), (pO|) and (|31|). The 
center of this algebra is infinite-dimensional and its structure depends on the 
dimension of x-space. In 2-dimensional case the center is formed by : 

Mm)) (54) 

see , for the discussion. In 3-dimentional case the general answer is not 
known to our knowledge and its construction is out of the scope of our paper. 
As an example we mention that Hopf invariant 

Q = Jd'xe,,,v,{x)VMS) (55) 

belongs to this center. 

4 Infinite-dimensional symmetry and inte- 
grals of motion. 

As we have mentioned above, the Lagrangian ([T^) possesses the invariance 
with respect to the "volume preserving" group of diffeomorphisms SDif f[D], 
where D O R^. Indeed, let us write the Lagrangian ([TtD in terms of the 
Lagrangian density 

L = [ d'^C{^,) (56) 



and consider the transformations from SDiff[D] of the coordinates S,i & D 

=0,(e.), (57) 
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det^-^ = 1. (58) 

Apparently, due to ( |55D we obtain: 

L = j d?<P{ii)C mi)) = I d'^C (59) 

and according to Noether's theorem this invariance results in the existence of 
an infinite set of integrals of motion. To obtain these integrals we first need 
to find the parametrization of the transformations (0), (|58D in the vicinity 
of identity transformation: 

U^^)=^j+a,i^,). (60) 
From (^) follows the equation for aj(^j): 



90 



0. (61) 



Further we must explicitly take into account that the volume preserving 
diffeomorphism (|57| ) leaves the boundary of D invariant. We shall limit 
ourself with the case when D is formed by extraction of the domain with the 
differentiable boundary given by 

gi^i) = (62) 

from R^. Physically that means that we put in the fiuid the fixed body, 
the shape of which is given by (|62D. The condition that the infinitesimal 
diffeomorphism (|60|) preserves D in this case is 

^?(e, + «,(e.))|^(^^)^„ = 0, (63) 



or 



«,(6)V,^?(e.) ^(^^^^^ = 0. (64) 

Geometrically equation ( |64D means that the vector a(^j) is tangent to the 
surface, defined by (§2]), because the vector Vjg{^i) ^ is proportional to 

the normal ^j(^i) of the surface (^) at the point ^j. 
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From Noether's theorem we obtain that the invariance of the Lagrangian 
with respect to the transformation (|60|) gives the following conservation law: 

1 / d'CPm{^^)^^^MC^)=0, (65) 

ot Jd o^i 

where ai{^i) satisfies to the conditions (|6lD and (|6^) . The existence of these 
conditions forbids to take the variation of the l.h.s. of ( |5BD over ai{^i) and 
obtain the local form of integrals of motion. For that we need to extract 
from (^) and ( |5iD the integral properties of az(^i). 

Consider an arbitrary, single- valued, differentiable in D function I3{^i). 
Then the following equations are valid: 



= / d'^mi) 

JD 



= - jj'^^M^f-^^ (66) 

The first equality is valid due to condition (|6l|). Using Stokes theorem we 
can transform the integral of total derivative in the last equality (]66|): 

Id ^'^a^(^(^^)«^(^^)) = L = 0' (67) 



where dD denotes the boundary of D. The last integral in (67) vanishes due 
to condition (0) because the differential dSj is proportional to the normal 
vector of the surface dD, defined by (^). From (|66[) and (|6^) we conclude 
that 

^'e«.te)^ = o (68) 

D Ol;j 

for any smooth, differentiable /?(^i). Taking this property of aj(^j) into ac- 
count we obtain from the conservation laws (^) that the quantities 

Jk{^^) = PU^^)^^^ (69) 

are conserved modulo some term which is the gradient of a scalar. In partic- 
ular, that means that 

= ^^M-^Jm = J- ipn.{i^f-^^] ■ (70) 
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is the integrals of motion. Note, that as the group of invariance is infinite- 
dimensional, we obtain an infinite number of integrals of motion. With re- 
spect to Poisson brackets ( PB| ) the i?j(^j)'s form an algebra. This algebra 
could be written in a compact form for integrated objects 

R[<P\ = I £^<P,i^,)R,{^,), (71) 

where are smooth, rapidly decreasing functions. The algebra of R[(f)] 

induced by Poisson brackets ( p3D has the form: 

R[ip]} = R[curl<p X curl^)] (72) 

The construction of the x-dependent object, corresponding to Rj{Ci) is not 
an easy task, because our "averaging" with S{x — x{$,i)) will introduce time 
dependence and instead of conserved object we shall obtain a density, whose 
time derivative gives a divergence of a "current". Therefore we need to in- 
troduce another kind of averaging without explicit refering to the x(^j) coor- 
dinates. For that recall that under diffeomorphism a closed loop transforms 
into closed loop. Then let us consider such a loop A and a surface a whose 
boundary is A. The integral 

V=jdS,R,{i.) (73) 

where the vector dSj is as usually the area element times the vector, per- 
pendicular to the surface, is conserved, because of the conservation of -Rj(^j). 
Further, from Stokes theorem we have: 

^ = //Ol'-(e.)^^. (74) 
Changing the variables in ( [71D we obtain: 

"^^^^'M-Za'^^"^'^-'- '''' 

where A is the image of the loop A under diffeomorphism The 
object (|75| ) is very well known in hydrodynamics as the "circulation" and 
its conservation is known as W.Thompson theorem The relation of the 
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circulation conservation with the invariance under special diffeomorphisms 
was first explicitly established in |Tl| , though it also could be extracted from 
general discussion in Appendix 2 of 

Conservation of circulation is not the only consequence of (|65|). Consider 
for example the case of 2-dimensional space. Here, instead of the conserved 
vector Rj{C,i) we shall have the conserved scalar 

This scalar defines the following integrals of motion: 

Changing variables in (|77D C,j — > C,j{xi) which is possible, because ^j{xi) is a 
diffeomorphism of D, we obtain: 

^dV(x.)(^eHAi(e(a;))Ar(e(a:))^^^^^)", (78) 

where the matrix Al.{C,) was defined in (|MD. We can present J„ as 

using the following property of 2-dimensional matrix AKC,): 

ekiAi{ax))AT{ax)) = e.mdetAi^x)) = e,m-^y (80) 

Recall that in the case of the Lagrangian (0), we are considering 

^^ = mv,ix.) (81) 

Therefore the integrals /„ coincide with the center (|5^) of the Poisson algebra 
331). 
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In the case of 3-dimensional space we can construct the analogous inte- 
grals of motion, integrating the products of the vector (ffOl) : 

%,i.....„ = (82) 

Changing variables — > Xi^ clS above we shall obtain: 

R,{0 - RAa^)) = ^Mk{m)m{^) f^f}''^^ - (83) 
In the 3-dimensional case the matrices A^{E,{x)) satisfy the equation: 

e,kiAt{ax))A^{i{x)) = emnr^^detAiax)), (84) 
Therefore (|83D takes the form: 

RM^)) = -±-e,.J-^?^ (85) 

p[Xi) dXr OXra 

The integrals Kn in (|82|) become 

= \^S'xp{x,)R,M^))R,M^))...R,Sm) (86) 

Apparently, due to the presence of in (|85|) , we can not in this case 

express only in terms of velocity and density, i.e. the Eulerian description 
does not admit this kind of integrals of motion. 

In the 3-dimensional case there is one more integral, which does not ex- 
ist in any other dimension. Recall that the vector Jki^i), given by (|69D is 
conserved modulo gradient, therefore the integral 

Q = f d^^M^.)Rk{^^) (87) 

is conserved because R{^i) = curlJ{^i). Transforming the ^j-dependent vari- 
ables to the Xj-dependent ones in (pT]) we obtain Hopf invariant: 

Q = d^e,uipM^)f-^^§^- (88) 
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5 Inclusion of the electromagnetic interac- 
tion. Plasma. 



We shall consider plasma as a fluid of two components, namely electrons with 
mass m and electric charge (— e) and ions with mass M and charge (+e). 
The coordinates of electrons we shall denote as x{C,i), while the coordinates 
of ions will be The interaction of the components of the plasma with 

the electromagnetic field is governed by the following Lagrangian: 



e[Ao{x{^i,t))~xi^,,t)A{xi^,,t)) 

-i J d?xF^\x,t)F^,{x,t), 
where Lq and Lq°" are "free" Lagrangians 



(89) 



T el 
^0 



•2 , 



- r{det- 



T ton 
^0 



MX{E,,t) 



- f^idet 



dXj{E„t) 



2 ' dEk 

and Ff^i^{x,t) denotes the electromagnetic field tensor: 
F^^{x,t) = d^A^{x,t) - d^A^{x,t) 



(90) 



(91) 



The advantage of Lagrangian description of fluid (plasma) is clear. Using 
of the coordinates of charged particles as the fundamental variables makes 
explicit the interaction with electromagnetic field. 

The Lagrangian ( ^9l) possesses U{1) gauge invariance. We shall use the 



usual Hamiltonian formalism for the constraint system [|T^. The canonical 
variables are: 



5L 

6x{^, 
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X{ 



A{xi 



(Xi 



6L 

6L 

6A{xi 
= 0, 



where E{xi) is the electric field strength: 

E{xi) = -VAoixi) - A{xi). 



(92) 



(93) 



The last of the equations (|92D is actually the primary constraint. The Legen- 
dre transformation of the Lagrangian (PW) gives us the canonical Hamiltonian: 



H = Jd^l {PLi^^) + H\x,)) + Ao(x,)Vn(a;, 



+ 



+ 



Jd'^MI (^(S.) - eA{X{E,))) + eAo{X{E,)) + /-(rfet^^; 



(94) 



where H{xi) = curlA{xi) is the magnetic field strength. The requirement of 
the conservation of the primary constraint Ilo{xi) = gives the secondary 
constraint ( Gauss law ): 



Vj {Pem{Xi))- - e {pel{Xi) - Pion{Xi)) = 0. 



(95) 



Now we can add to the primary constraint no(xj) = the gauge fixing condi- 
tion Ao{xi) = and eliminate these variables from consideration. Introducing 
the x-dependent functions instead of u and [/-dependent, as we did in the 
2-nd section we obtain the Hamiltonian of the plasma in the following form: 

(l{xi) + epeiixi)A{xi)] 



(^L{xi) + epion{xi)A{xi) 



2mpei{xi) 



2mpion{Xi) 

where we have introduced 

l{Xi) 

Lixi) = 



+ Vel{Pel{Xi)) + Vion{pion{Xi)) 



= id'^msix-£m 

Jd'EP{E,)6{x-X{E,)) 



(96) 



(97) 
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This Hamiltonian is gauge invariant with respect to the transformations, 
generated by the constraint (|95D . Further we impose the Coulomb gauge 
condition on the electromagnetic field 



V,A,{xi) = 



(98) 



and following the usual procedure will eliminate the longitudinal components 
of A{xi) and Pem{xi). As a result of the gauge fixing, the longitudinal part 
of Pem{xi) give rise to the Coulomb term in the Hamiltonian 



Col 



(fx {pel{Xi) - Pion{Xi)) 



4:7r\x — y\ 



{peliVi 



and the whole Hamiltonian takes the following form: 

H = Hcoi + Id^x[l (P2^, (x,) + HHxt 
(l{xi) + epei{xi)A±{xi) 



PioniVi)) ■ 

(99) 



, 2 



2mpei{x,i) 
(^L{xi) + epion{xi)A_i{x, 



2mpi 



[Xi 



+ Vel{pel{,Xi)) 



~\~ '^ioni^Pioni^Xi^) 



(100) 



where the subscript ± denotes the transverse components of electromagnetic 
variables, for which the Poisson (Dirac) brackets are given by |jl2| 



{PlJ{x),A^,{y)} 



1 

A 



djdk ]6{x-y). 



(101) 



As usually, the Coulomb term (^11) contains an infinite additive part of 
"self interaction" which should be subtracted. In equilibrium plasma the 
Coulomb interaction is known to be screened by the cloud and the residual 
Debye interaction is a short rang one (see e.g. ll^). In Jl^ it was stated that 
Coulomb term is reduced to the local functional of the difference of charge 
densities {pei{xi) - Pion{xi)). 

Following our consideration of a fiuid in the Section 3, we can introduce 
the canonical coordinates ^{xi),7f{xi) for electron and S(xj),H(xj) for ion 
components of the plasma: 



{T^j{xi),^kiyi)} = 5jkS{x-y), 
{Uj{xi),Ek{y^)} = 6jk6{x-y), 



(102) 
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for which 



L,(x,)=S(x,)^^. (103) 

Substituting ( |103| ) into ( |100| ) we obtain the Hamiltonian in terms of the 
canonical variables. 

The Lagrangian (^9]) is invarian with respect to the volume preserving 
diffeomorphisms of both components of plasma separately, so we shall have 
in this case two sets of conservation laws - one for electrons, the other for ions. 
Applying the procedure, which we have described in the previous section we 
shall construct the conserved circulations: 

n-' = kdx,^, (104) 

where lj{x) and Lj{x) are given by ( |103| ). In the same way we can construct 
the analogues of the integrals (^61) and (|88|) for this case. Note that in the 
case of plasma the equation (|SB| ) is not valid due to the presence of the 
electromagnetic field, instead we have: 



Kx) = peiix) [mvei{x) - eA{x) ^ 

) (Mvionix) + eA{x)) , (105) 

therefore the equations (|104|) are not the circulations of the velocities. In the 
same time, adding circulations and V^""" on the common contour A we 
obtain: 

Vk = Vt + = £ dxj (mveiix) + Mvionix)) , (106) 

so this linear combination of circulations of electons and ions velocities is 
conserved even in the presence of the electromagnatic interaction. 

Summarizing we can say that the phase space of plasma F in Coulomb 
gauge is the space with coordinates A±{x), Pl^i^x); ^{xi),'7i{xi);E{xi),Il{xi) 
with Poisson brackets given by (|101|) and (|102|) . The evolution of a state in 
r is difined by the Hamiltonian ( |100| ). 
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6 Concluding remarks 



The Hamiltonian approach for ideal fluid and plasma, considered in the 
present paper in some aspects resembles the Hamiltonian formulation of a 
classical field. It deals with the system with infinite number of degrees of 
freedom - "particles" which constitutes a continuous media. These particles 
interact with nearest neighboors due to the potential part of the Lagrangian 
or Hamiltonian, as the field amplitudes in conventional field theory. The dif- 
ference with the later, which we emphasized in the course of the paper is that 
the constituents of the fluid change their position in space while evolution 
and to have a local description of fluid in terms of x-dependent variables we 
need the projection realized by e.g. equation (^). The fluid dynamics could 
be also compared with the mechanics of the extended objects like n-branes. 
These objects are also n-dimensional media embeded in the space-time, but 
its constituents are indistinguishable and due to this property the general dif- 
fiomorphisms of its coordinates appear as the gauge invariance of the theory. 
In the case of the fluid it is supposed that the constituents are distinguishable 
and, as a result only special diffiomorphisms preserve the Lagrangian giving 
rise to the Noether's integrals of motion. 

So far we have considered only classical mechanics of the fluid and plasma, 
but as we know any classical Hamiltonian system could be quantized. Here 
as a quantization we mean a formal procedure of the transition from from 
classical canonical coordinates to the operators, satisfying commutation re- 
lations inherited from Poisson brackets and the construction of a represen- 
tation of these operators in an appropriate space. Of course this procedure 
has no physical meaning in the case of plasma, but for a fluid it could be 
quite reasonable providing us with the theory of an ideal quantum fluid. The 
quantum fluid could be considered as the new representation of a quantum 
field in which the eigenvalue of the operator of density of particles may con- 
tinuos function of x. The matter is that in the case of usual quantum field 
we know only Fock representation |]T3|, for which the operator of density of 
particles in space (or in momentum space) has the eigenvalues which are 
the superposition of 5-functions, that means that the average density is zero 
and to consider the states with finite density we need to go outside the Fock 
space. 
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